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Abstract 

We  develop  a  computer  program  for  Kalman  filtering,  smoothing,  extrapolation 
in  dose-response  experiments,  based  on  the  method  proposed  by  Meinhold  and 
Singpurwalla  in  1987,  and  we  apply  this  program  to  some  published  data  on  doses  of 
bischloromethyl  ether  administered  to  rats,  discussed  in  Chen  and  Singpurwalla 
(1989). 
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1 .  Introduction 


In  a  biological  or  an  engineering  system,  we  are  often  required  to  infer  the 
response  at  a  very  low  dose  (or  stress)  level,  based  on  the  data  of  the  observed 
values  of  responses  at  higher  dose  (or  stress)  levels.  Meinhold  and  Singpurwalla 
(1987)  have  proposed  a  method  involving  a  use  of  the  technology  of  Kalman  filtering 
and  Kalman  filter  smoothing,  for  inference  and  extrapolations  in  certain  dose- 
response,  damage-assessment  and  accelerated  -  life-testing  studies.  Singpurwalla  and 
Chen  (1989)  have  discussed  issues  pertaining  to  a  practical  implementation  of  this 
methodology.  Aiming  at  an  easy  application  of  the  methodology  to  real  cases,  we 
have  developed  a  computer  program,  which  is  described  herein. 


2.  Algorithm 

2.1  Model  (Meinhold  and  Singpurwalla  (1987)) 

Assume  that  E(Y(x))  =  exp  (-  a(x)  x  where  a(x),  /3(x)  >  0  and  ‘•E(Y(x))' 

denotes  “the  expectation  of  response  Y  at  dose  (or  stress)  level  x  .  Suppose  that 
Y*(x)  log(-log  Y(x))  and  Y*(x)  ^  .N-(Mx),  ^^(x)).  Then,  the  observation  equation 

of  the  Kalman  filter  is 

Y-  -  Fj  (2-1 

the  system  equation  of  the  Kalman  filter  is 


where  Y*  =  Y*(Xj),  0j  =  d  log  Xj]  = 

i/ ■  ~  jV(0,  V  )  and  W;  —  JV(0,  W  J;  u :  is  assumed  independent  of  Wj;  Xj_^ 
J  J  J  j  J 

(stress)  previous  to  dose  Xj. 


is  the  dose 


Computational  Method  (Singpurwalla  and  Chen  (1989) 


Assume  that  ~  ^ j_i).  then,  upon  observing  Y*,  0^  ~  -Y(0j,  tj), 


where 


«j  -  »j-i " 


“  ^j-1 


(2.3) 


X  (I  - 


where  I  is  an  identity  matrix.  The  relationship  in  (2.3)  is  referred  to  as  the 
forward  recurrence  equations  of  the  Kalman  filter.  The  backward  recurrence 
equations  are 

§j(T)  -  9j  +  Jj+i  IVl<T)  -  ¥ 

h/T)  -  Ij- 

where  J  =  S  i  R'\  0  (T)  is  an  estimation  of  0:  based  on  all  data,  and  we  assume 
j  j-l  J  J 

that  0-  ~  WjlT),  S:(T))  during  backward  recurrence. 

J  J  J 

Once  0-j-(T)  is  obtained,  inference  about  Y(xj^P  ^^T+1  ®  (stress) 

lower  than  x^-)  is 

E(Y(xj^p)  ^  exp(-  exp  0j(T))),  (2.5) 


since  Y*(x-p^P  ^  0'j'(T),  T'j'(T)  F'p^p. 

Besides,  to  implement  the  Kalman  filter,  we  are  required  to  give  some 
necessary  starting  values,  say  starting  doses  (stresses)  x(0)  and  x(-l),  assumed 


responses  y(x(0))  and  y(x(-l))  for  getting  assume  variances  of  y*(x(i)) 

=  log  (-log  y(x(i))),  i=0,  -1,  and  their  correlation  for  calculate  Jq.  We  also  propose 
W(Xj)  =  .5  ♦  (Xj_j  -  Xj)^  •  SCxq).  (2 

V(Xj)  =  (1/exp  (exp(Fj0j_p)r(l  -  Y  exp(exp(Fj0j_2)))"'.  (2 


3.  Usage  of  the  Computer  Program 

The  program  has  been  written  in  Mikrosoft  QuickBASIC  version  4.0,  to  be 
run  on  an  IBM  Personal  computer  or  “compatible”.  A  printer  should  be  available. 

The  filename  of  the  program  is  tr896.BAS. 

This  program’s  flow  chart  is  drawn  up  below. 


END 

Figure  3.1  Flow  Chart 
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This  program  contains  a  subroutine  named  “ENTRY”  which  allows  either  to 
input  data  and  create  a  data  file  or  to  retrieve  input  from  an  existing  file.  As  soon 
as  the  program  gets  started,  "ENTRY”  will  be  called.  At  this  moment  figure  3.2 
appears  on  the  screen. 

FILTERING,  SMOOTHING,  .4ND  EXTRAPOLATIONS 
IN  DOSE-RESPONSE  EXPERIMENTS 

An  Input  File  is  required  for  the  calculations 

PRESS  R  to  retrieve  input  from  existing  file, 

C  to  create  new  input  file,  or 

Q  to  quit. 

Figure  3.2  Selecting  data-input  form 

There  are  three  options; 

1)  to  quit  and  return  to  the  system,  press  key  “Q”; 

2)  to  input  data  interactively  and  generate  a  data  file,  press  key  “C  ; 

3)  to  retrieve  input  from  an  existing  file,  press  key  “R”. 
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screen 


Once  key  "C”  is  pressed,  the  user  should  input  the  data  requested  on  the 
(see  figure  3.3).  The  input  data  will  be  saved  in  an  ASCII  data  file. 


Total  number  of  pairs  (stress  level  -  response) 
considered  by  the  experiment  (N'D)  . 

Total  number  of  pairs  (stress  level  -  response)  for  which 
predictions  are  required  (NP)  . . . 

Starting  value  of  the  stress  level  (x(0)) . 

Stress  level  Observed  response 

1.  - 

2.  - 

“ND.”  - 
Stress  levels 

1.  - 

2.  -  ? 


? 

9 


? 

9 


?  ? 
for  which  responses  are  to  be  predicted 


“ND.”  -  ? 

Second  starting  value  of  the  stress  level  (xC-O) . 

Assumed  response  at  the  starting  stress  level  x(0)  . 

Assumed  response  at  the  second  starting  stress 

level  x(-l)  . 

Let  YO  correspond  to  Y  star  of  X  subscript  0  on  the  paper 

Assumed  variance  of  YO  . * . 

Let  Y1  correspond  to  Y  star  of  X  subscript  -1  on  the  paper 

Assumed  variance  of  Y1  . 

Assumed  correlation  of  YO  and  \  1  . 

File  name  where  input  data  will  be  saved  . 


Figure  3.3  Hints  for  data  input 
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Once  key  '‘R”  is  pressed,  the  user  should  enter  the  name  of  the  data  file  to 
be  retrieved,  which  format  should  meet  the  requirement  of  the  program,  i.e.  the 
sequence  of  the  data  should  be: 

the  number  of  data  paris  of  observed  responses  and 
correspondent  dose  (stress)  levels  (ND), 

the  number  of  dose  (stress)  levels  at  which  prediction 
is  required  (NP), 

starting  dose  (stress)  level  (x(0)), 

dose  (stress)  level  l(x(l))  and  observed  response  at 
level  x(l)  (Y(x(ll)), 

dose  (stress)  level  2  (x(2))  and  correspondent  response 
(Y(x(2))), 

dose  (stress)  level  ND  (x(ND))  and  correspondent  response 
(Y(X(ND))), 

X(ND  +  1)  -  the  1st  level  at  which  prediction  is  required, 

X(ND+NP)  -  the  last  level  at  which  prediction  is  required, 
another  starting  dose  (stress)  level  (X(-l))  (note  that 
X(i)  <  X(i-l),  i=0,  ,  .  .  ,  ND  +  NP,  is  required), 

assumed  responses  at  levels  x(0)  and  x(-l)  (Y(x(0)),  Y(x(-1))), 

assumed  variances  of  £n(-£n  Y(x(0)))  and  £n(-£nY(x(-l))), 

assumed  correlation  of  £n(-£nY(x(0)))  and  £n(-£nY(X(-l))). 

After  the  user  has  either  input  all  the  data  or  retrieved  the  data  from  an 
existing  file,  the  program  prints  all  input  data  (as  shown  in  Section  4,  via  an 
example)  and  starts  computing.  Finally,  computation  results  are  printed  out.  The 
print  out  is  composed  of  three  parts  -  estimation  of  parameters,  dose-responses 
(including  observed,  predicted,  filtered  and  smoothed  responses)  and  prediction  of 
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responses  at 


low  dose  levels  (as  shown  in  Section  4,  via  the  example. 


4.  An  Application  of  the  Program  to  Data  on  Respiratory  Tumor  of  Rats 


4.1  A  Croup  of  Data  on  Respiratory  Tumor  of  Rats 
and  a  Set-up  of  Starting  Values 

The  data  on  respiratory  tumor  of  rats  are  abstracted  from  the  Final  Report 
of  the  SCFSC  (1980),  which  are  presented  in  Table  3.1. 


Table  3.1 

Dose-Response  Data  on  Respiratory  Tumor  of  Rats 


Dose:  x^ 

10  20  40  60  80  100 

Observed 

Response 

40/41  43/46  14/18  14/18  19/34  8/20 

In  the  table,  datum  “40/41”  denote  the  fact  that  40  out  of  41  rats  have  not 
developed  tumors. 

Our  starting  values  are  given  as  following: 

x(0)  =  250,  x(-l)  =  280,  y(x(0))  =  .005,  y(x(-l))  =  .001, 

Var  (y*(x(0)))  =  .OOl^,  Var  (y*(x(-l)))  =  .003^,  p(y*(x(0)),  y'(x(-l)))  =  . 

The  doses  at  which  response  prediction  is  wanted,  are  7,  5,  2,  1. 


/ 


4.2  The  Printout  of  Input  Data  and  Computation  Results 


FILTERING,  SMOOTHING,  AND  EXTRAPOLATIONS 
IN  DOSE-RESPONSE  EXPERIMENTS 

INPUT  DATA  SAVED  AT  “FILENAMES”: 

Total  number  of  pairs  (stress  level  -  response)  considered 
by  the  experiment  (ND) 

Total  number  of  pairs  (stress  level  -  response)  for  which 
predictions  are  required  (NP) 

Starting  value  of  the  stress  level  (X(0)) 

Second  starting  value  of  the  stress  level  (X(-l)) 

Assumed  response  at  the  starting  stress  level  X(0) 

Assumed  response  at  the  second  stress  level  X(-l) 

Let  YO  correspond  to  Y  star  of  X  subscript  0  on  the  paper. 
Assumed  variance  of  ln(-ln(Y0)) 

Let  Y1  correspond  to  Y  star  of  X  subscript  -1  on  the  paper. 
Assumed  variance  of  ln(-ln(Yl)) 

Assumed  correlation  of  ln(-ln(\0))  and  ln(-ln(Yl)) 


Stress  levels 

Observed  responses 

100 

.4 

80 

.5588 

60 

.7778 

40 

.7778 

20 

.9348 

10 

.9756 

Stress  levels  for  which  responses  are  to  be  predicted: 

7 

5 

2 

1 


6 

4 

250 

280 

1 ,667389 
1 .932645 

.000001 

9E-08 

2.1E-07 
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From  the  starting  conditions  the  following  are  the  parameters 

ALPHA,  GAMMA,  BETA  and  the  covariance  matrix  of  GAMMA  and  BETA, 
all  of  them  at  the  starting  stress  level  X(0): 

START.ALF=  1.292872E-05  START.GAM=-1 1.25606  START.BET=  2.340585 
SIGMAO:  .0026858  8.486889E-05  -4.774242E-04 


FILTERING,  SMOOTHING,  AND  EXTRAPOLATIONS 
IN  DOSE-RESPONSE  EXPERIMENTS 
PROGRAM  RESULTS 

TABLE  1.  Estimation  of  Parameters 


DOSE 

FILTERING 

SMOOTHING 

GAMMA 

BETA 

GAMMA 

BETA 

100 

-11.24675 

2.33893 

-11.22904 

2.335781 

80 

-11.21696 

2.333635 

-11.17025 

2.325325 

60 

-11.14165 

2.32025 

-10.1507 

2.144237 

40 

-9.543344 

2.036153 

-8.848135 

1.912623 

20 

-8.113652 

1.782019 

-8.044925 

1.769795 

10 

-7.605752 

1.691735 

-7.605752 

1.691735 

TABLE  2.  Dose 

-  Responses 

DOSE 

OBSERVED 

PREDICTED 

FILTERED 

SMOOTHED 

100 

.4 

.5376881 

.5371238 

.5360512 

80 

.5588 

.6915619 

.6898784 

.6872347 

60 

.7778 

.8271979 

.8239535 

.7758961 

40 

.7778 

.9272031 

.877179 

.8466142 

20 

.9348 

.9685547 

.9395663 

.9376571 

10 

.9756 

.9820372 

.9758289 

.9758289 

TABLE  3.  Prediction  of  Responses  at  Low  Levels 


DOSE 

PREDICTION 

7 

.986706 

5 

.992454 

2 

.998393 

1 

.999502 
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